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1  Introduction 


Many  facets  are  involved  in  the  design  of  a  spread  spectrum  communica¬ 
tions  system  [23].  One  important  consideration  is  determining  a  method 
to  best  recover  the  transmitted  signal  when  it  is  subjected  to  interference 
in  the  transmission  path.  This  interference  is  often  highly  correlated  and 
not  necessarily  Gaussian  as  in  typical  interference  models  [Ij.  The  linear 
correlator  realization  of  the  matched  filter  is  no  longer  optimal  for  this  in¬ 
terference  environment.  Locally  optimum  (LO)  detection  provides  a  method 
for  circumventing  this  problem  if  the  probability  density  function  (pdf)  of 
the  interference  is  known.  However,  in  cases  where  the  interference  exhibits 
strong  self-correlation,  traditional  LO  methods  exhibit  poor  performance  [3]. 

LO  detectors  with  memory  more  successfully  combat  this  type  of  distur¬ 
bance.  The  disadvantage  of  a  memory-based  processor  is  the  rapid  increase  in 
the  dimensionality  of  the  joint  probability  density  function  (jpdf)  noise  vector 
making  the  LO  detector  nonlinearity  unwieldy  [2].  To  alleviate  this  problem, 
frequency  domain  methods  are  used  to  determine  a  -order  autoregressive 
(AR(P))  model  of  the  channel  disturbance.  The  AR  methodology  reduces 
the  dimensionality  of  the  underlying  jpdf  to  that  of  the  model  order.  The 
statistics  of  the  jpdf  remain  an  issue.  To  this  end,  probability  density  function 
(pdf)  estimation  techniques  for  independent  identically  distributed  (iid)  noise 
samples  prove  useful  (11)  since  the  input  sequence  to  the  AR  model  is  iid, 
albeit  characterized  by  some  unknown  pdf. 

The  necessary  AR  model  parameters  must  be  estimated  by  well-known 
spectral  estimation  techniques  [1].  Considering  that  each  received  symbol 
will  be  limited  by  the  sampling  rate,  an  estimation  technique  must  provide 
accurate  results  for  a  finite  number  of  observations  of  the  random  process. 
The  Yule- Walker  Method  [1]  works  best  when  the  number  of  samples  ap¬ 
proaches  infinity.  Other  techniques  [4]  can  estimate  the  model  parameters 
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given  very  few  samples. 

This  report  presents  the  derivation  and  simulation  of  a  robust  locally  op¬ 
timum  (LO)  nonlinear  spread  spectrum  receiver.  The  signaling  environment 
consists  of  the  desired  received  signal  in  correlated  interference  and  thermal 
noise.  Autoregressive  (AR)  spectral  modeling  methods  and  a  histogram  ap¬ 
proximation  of  the  probability  density  function  are  employed.  The  system 
under  consideration  is  a  binary  phase-shift  keyed  (BPSK)  communications 
system.  Results  from  simulations  [9]  indicate  that  applying  AR  modeling 
to  the  environment  improves  performance  substantially  in  most  cases.  This 
document  summarizes  the  work  to  date. 


2  The  Locally  Optimum  Detector 

Consider  the  case  of  binary  detection  of  real,  discrete  signals  in  additive 
noise.  In  binary  signal  detection  the  goal  is  to  decide  which  of  two  possible 
information  signals  is  present  at  the  detector.  For  example,  in  a  binary 
phased  shift  keyed  (BPSK)  communications  system  [22],  the  receiver  must 
decide  whether  a  value  of  -fl  or  —1  was  sent  by  the  transmitter.  Stated 
more  formally,  the  goal  of  the  detector  is  to  correctly  choose  one  of  the 
following  two  hypotheses: 

Hi  :  Signal  Si  present  (1) 

Hq  :  Signal  Sq  present, 

where  Si  and  Sq  are  vectors  and  the  notation  x  =  [xi---xnY  denotes  a 
column  vector  of  length  N.  The  value  of  N  is  the  number  of  signal  samples 
in  a  given  observation  period. 

The  derivation  of  the  corresponding  LO  detector  with  memory  is  sum¬ 
marized  by  Jacklin  and  Ucci  in  Chapter  9  of  [11],  and  is  repeated  here  for 
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convenience.  To  begin,  let  the  received  signal  vector  of  length  N  be  given  by 


r  —  Sm  "b  H ) 


where  m  =  0  or  1  and  n  is  a  random  noise  vector.  Given  the  observation 
T  =  p,  the  optimum  detector  is  of  the  form  [17] 


choose  Hi 


ln[A(p)] 


choose  Hr 


where  the  log  likelihood  ratio  is  given  by  [17] 


ln[A(p)]  =  In 


Mp  I  Hi) 

Up  I  Ho) 


/n(P  -  Si) 

/n(p  -  So) 


assuming  that  s^  and  n  are  independent.  Approximating  Eq.  (4)  using  a 
first-order  Taylor  series  expansion  yields  the  following  result: 


^  o 

InA(p)]  «  “  ■soi)  ^ln[/„(p)] . 

i=l 

^.From  Eq.  (5),  the  resulting  decision  statistic  becomes 


where 


^(P)=  “  Soi)9i{p), 


^  ap.  “  Up) 


is  called  the  LO  nonlinearity  with  memory.  Thus,  the  LO  detector  with 
memory  for  the  binary  detection  of  known  signals  in  additive  noise  is 

choose  Hi 

Kp)  =  -  ■Soi)5i(p)  7.  (8) 


choose  Hr 
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with  gi{p)  given  by  Eq.  (7). 

If  n  is  a  first-order  Markov  random  sequence  [24],  then 

fm {m  I  ^i-2,  •  •  • ,  >7i)  =  /ni [Vi  I  Tli-i)  ■  (9) 

Using  Eq.  (9),  the  joint  pdf  of  n  can  be  reduced  to 

N 

fniv)  ~  n  fnii'Oi  I  Vi-l)  >  (10) 

<=1 

where  fmivi  I  ^o)=/ni(^i)-  The  LO  nonlinearity  requires  /n(p),  the  noise 
pdf  evaluated  at  the  value  of  the  received  signal.  From  Eq.  (10)  this  is  seen 
to  be 

fn{p)  =  lifmipilpi-l)-  (11) 

1=1 

Substituting  Eq.  (11)  into  Eq.  (7),  the  LO  nonlinearity  with  memory,  gi{p)  — 
-^ln[/„(p)],  becomes 


9iip)  =  ( 


[  “4: I  +Hfni+iiPi+i  I  pi)]) , 

for  i  =  1, ...,  A'  —  1 


(12) 


I  “^M/n/v(PiV  I  PiV-l)], 


for  i  =  N . 


Thus,  for  first-order  Markov  noise,  the  original  AT- variate  function  pi(p)  re¬ 
duces  to  a  function  of  at  most  three  variables:  pj-i,  pi,  and  pi+i. 


3  The  LO  Detector  for  Auto-Regressive  Noise 

The  focus  of  this  effort  is  the  development  of  robust  LO  detection  algorithms 
for  general  AR  noise  environments.  Suppose  the  noise  sample,  rij,  at  a  dis¬ 
crete  time  instant,  i,  can  be  modeled  as  a  -order  Markov  process.  Then 
the  autoregressive  noise  model  is  [3] 
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(13) 


0, 


for  z  e  [  1,  N], 
for  z  <  0 


where  the  {aj  are  the  AR(P)  coefficients  and  Wi  represents  an  indepen¬ 
dent  identically  distributed  (iid)  random  process  and  the  vectors  contain  N 
samples.  In  this  case,  the  noise  pdf  is, 


N 


fniv)  =  fniVu  ■,Vn)  =  '[[  fnM  I  Vi-U 


•  •  •  1  Vi-p) 


(14) 


where, 


Imirn  1  Vi-ly  •  •  •  )  Vi-p)  =  friiiVi)  (15) 

friiiVi  I  Vi-1,  1  Vi-p)  =  fuiiVi  I  Vi-1,  ,Vi)  for  z  =  2,  •  •  •  , P 

If  a  ’’block  approximation”  [1]  ^  is  assumed,  then  fn{v)  can  be  written 


N  !  p 

faiv)  ~  IT  ~  ^jVi-j 

i=l  \  j=0 

where  is  the  pdf  of  the  white  noise  process  and  oo  =  -1-  For  this  case, 

the  form  of  the  nonlinearity  is. 


min  (P.N  — i) 

E 

/=o 


'/»(»')  =  E  I  “  E  I  for  z  e  [  1,  jV ],  (17) 

j=0 


where, 


(18) 


0,  for  A:  <  0, 

Pk  foTke[l,N], 
and  h'{(jj)  is  the  derivative  of  h{u)),  and 

iln  general,  m  depends  on  the  last  “F”  previous  samples,  which  are  assumed  to  be 
zero  so  that,  771  =  wi.  This  only  affects  a  small  fraction  of  the  terms  for  iV  »  P. 
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/i(a))  =  ln[/^(£j)]  (19) 

Now  we  have  reduced  the  nonlinearity  to  Eq.  (17)  providing  an  intuitive 
understanding  to  the  ensuing  calculation.  Eq.  (17)  is  also  quite  suitable  for 
a  vector  representation  since  the  argument  of  h'(x)  can  be  represented  as 
a  dot  product  between  the  AR  parameters  and  a  received  vector.  In  turn, 
h'{x)  is  itself  a  matrix  of  these  dot  products  multiplied  by  the  vector  of  AR 
coefficients  thus  enabling  the  computation  of  the  necessary  nonlinearity  for 
the  symbol  of  interest. 

4  Histogram  and  AR  Estimation  Techniques 

Essential  to  the  computation  of  the  LO  detector  nonlinearity  are  the  noise 
pdf,  U  (w),  and  the  autoregressive  coefficients,  {oi}.  Since  neither  the  source 
statistics  of  the  driving  white  noise,  nor  the  AR  coefficients  are  known  a  pri¬ 
ori,  they  must  be  estimated. 

The  estimation  of  the  pdf  is  performed  via  a  histogram  method  employed 
in  the  LO  detector  without  memory  [1].  Pure  noise  data  is  not  available,  so 
the  received  data  must  be  used.  Then,  a  three-point  derivative  is  used  to 
compute  h'{p).  The  histogram  approach  is  chosen  because  of  its  ease  and 
simplicity  of  implementation  and  generally  acceptable  performance. 

The  AR  coefficients  are  estimated  using  the  modified  covariance  algorithm 
(MCA)  [12]  as  developed  [4]  and  simulated  and  tested  earlier  [5].  Recall  that 
the  MCA  minimizes  the  sum  of  the  forward  and  backward  prediction  errors: 

(20) 

where  the  forward  prediction  error  is: 
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(21) 


P 

(n)  =  x{n)  +  ^  ap{k)x{n  -  k) 

k=l 

and  likewise  the  backward  prediction  error  is: 

e~ {n)  =  x{n  - p)  + a*p{k)x{n  -  p  +  k)  (22) 

jfe=i 

and  the  error  is  given  by: 

=  ElKWP  +  (23) 

n=:p 

This  minimization  of  error  leads  to  the  normal  equations  that  are  central 
to  the  algorithm. 

^[rx{l,k)+rx{p-k,p-l)]ap{k)  =  -[rxil,0)+rxip,p-l)];l  =  (24) 

k=l 

where: 

ra;(/,  A:)  =  ^x(n  -  A:)a;*(n  -  0  (25) 

n—p 

and  the  {op}  vector  gives  the  estimates  of  the  AR  model  parameters  as 
required.  This  is  a  non-windowing  method  similar  to  the  covariance  method 
[12],  but  it  differs  from  the  latter  in  that  it  minimizes  the  sum  of  the  squares 
of  the  forward  and  backward  predictor  errors.  The  MCA  is  chosen  since, 
when  compared  to  other  AR  estimation  methods,  it  often  provides: 

•  stable,  high  resolution  spectral  estimates  for  a  small  number  of  samples 

•  exhibits  lower  sensitivity  to  phase  and  decreased  peak  shifting,  and 

•  is  not  subject  to  spectral  line  splitting. 
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The  Burg  Method  gives  the  benefit  of  providing  the  calculation  without 
requiring  the  estimation  of  the  autocorrelation  of  the  sequence  of  interest 
and  can  give  results  similar  to  the  Modified  Covariance  Algorithm.  The  Burg 
Method  would  be  preferable  in  an  embedded  system  with  a  finite  amount  of 
memory,  such  as  a  digital  signal  processor. 

5  System  Code  Translation 

To  make  the  computation  of  the  nonlinearity  viable  and  efiicient  in  an  embed¬ 
ded  system,  a  translation  to  C  or  C-b- 1-  is  an  essential  step  toward  real-time 
integration  on  a  digital  signal  processor  (DSP).  Two  possible  implementa¬ 
tions  were  pursued:  one  using  ANSI  C  and  the  other  using  C-I-+  [7].  In 
the  translation  to  C,  the  core  modules  of  filtering  and  histogram  generation 
were  successfully  completed.  The  C-f-+  translation  presents  a  simulation  of 
a  typical  DSSS  communication  system  for  the  variety  of  detectors  considered 
in  Section  6. 

Currently,  the  simulation  system  exists  only  in  MATLAB©  format.  To 
port  it  to  C/C-b+  a  translation  of  individual  MATLAB  modules  is  neces¬ 
sary.  The  first  module  translated  was  based  on  MATLAB’s  filter ()  function 
which  implements  an  infinite  impulse  response  (HR)  filter  [22].  We  imple¬ 
mented  the  direct  form  II  structure  because  of  its  simplicity,  minimal  memory 
requirements  and  fast  performance.  This  structure  requires  M  +  N  +  1  mul¬ 
tiplications  and  the  number  of  delays  is  of  the  order,  N,  where  M  and  N  are 
the  number  of  zeros  and  poles,  respectively,  of  the  filter  of  interest. 

The  filter 0  function  takes  three  null-terminated  arrays  as  inputs:  b[  ], 
representing  the  numerator  coefficients,  a[  ],  representing  the  denominator 
coefficients,  and  datain[  ]  representing  the  input  data  vector.  The  output  is 
a  single  vector,  y[]. 
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The  next  translated  function  was  the  MATLAB  hist()  function.  For  our 
purposes,  the  histogram  function  must  perform  the  following  tasks:  return 
the  distribution  and  bin  centers  of  an  input  array  among  a  specified  number 
of  bins. 

Once  implemented,  both  functions  were  tested  and  compared  to  results 
from  a  M.ATLAB  simulation.  In  testing  filter. c,  a  wide  range  of  varying 
inputs  and  filter  coefficients  were  examined.  Integer  and  fractional  input 
values  (positive  and  negative)  were  tested  as  well  as  varying  vector  lengths. 
A  filter  that  was  implemented  along  with  its  input  sequence  are  given  as: 

•  6[]  =  {1.22, -2, -2.01, -4.87} 

•  a[]  =  (4, 2, 3,  -1.54, 5.98,  -0.55} 

•  dataJn[  ]  =  (3, 9.2, 5.11, 6,  -6, 12, 0.77, 0.95,  -3.3344} 

Using  MATLAB’s  filter ()  function  as  a  baseline,  the  results  of  passing  the 
above  input  signal  through  filter,  c  revealed  virtually  perfect  agreement  [7]. 

In  testing  histogram,  c,  a  wide  range  of  bin  numbers  and  data  values  were 
examined.  Similar  to  the  filter,  c  case,  fractional  and  integer  data  were  used. 
To  illustrate  the  performance  of  the  C  implementation,  the  following  input 
signal  was  passed  through  the  histogram  function  which  was  divided  into 
binJnfo  =  5  bins: 

signall  ]  =  (10000, 85, 0, 0, 5,  -66.7, 88, 90, 0.555, 3.66,  -9.999, 4.665,  -25.6} 

When  compared  with  the  results  of  the  MATLAB  simulation,  again,  vir¬ 
tually  perfect  agreement  was  found.  The  results  were: 

6m_center.s  =  {939.97,2953.31,4966.65,6979.99,8993.33} 

f  requency -Count  =  {12,0,0.0,1} 
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The  C++  simulation  implemented  the  full  DSSS  communication  system. 
C++  routines  were  written  to  parallel  the  MATLAB  simulation  for  the  de¬ 
tectors  analyzed  in  Section  6.  An  object-oriented  approach  to  system  design 
was  utilized.  This  allows  for  a  higher  degree  of  flexibility  [21].  The  class 
hierarchies  for  the  data  structures  were  developed  [8]. 

The  base  class  for  all  of  the  signals  is  Matrix.  The  class,  Signal,  is  a 
special  case  of  the  Matrix  class  with  some  additions.  The  Signal  class  consists 
of  the  BPSK  Signal,  the  Noise  Signal  and  the  Interferer  Signal.  The  Noise 
Signal  represents  the  background  thermal  noise  process  and  the  Interferer 
Signal  represents  one  of  a  library  of  different  types  of  high-power  interferers. 
Currently,  the  only  two  interference  types  are  partial  band  and  continuous 
wave. 

The  translation  used  Template  classes  allowing  the  simulation  to  be  per¬ 
formed  in  both  ’’float”  and  ’’double”  data  types.  This  allows  an  examination 
to  be  made  of  the  trade-off  between  precision  and  efficiency. 

Simulations  were  executed  for  the  linear  correlator  for  ISR’s  of  2dB,  5dB, 
lOdB,  20dB,  and  25dB.  The  error  performance  observed  was  as  expected 
i.e.  as  interferer  power  decreased,  so  did  the  probability  of  bit  error.  This 
paralleled  the  MATLAB  simulation  results.  Run  times  of  the  entire  C++ 
simulation  were  much  faster  than  the  MATLAB  code.  For  a  partial  band 
interferer,  the  relative  performances  of  the  five  detectors  also  matched  the 
MATLAB  simulation  results. 

6  Simulation  Results 

The  following  figures  indicate  the  performance  of  the  ARLO  detector  versus 
four  other  detectors:  a  linear  correlator,  the  LO  detector,  a  linear  correlator 
and  an  LO  detector  both  preceded  by  an  AR  whitening  filter.  The  inter- 
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ference  environment  consists  of  of  any  one  of  the  following;  a  partial  band 
interferer,  a  CW  tone,  multiple  CW  tones  or  a  mixture  of  partial  band  and 
CW  interference.  The  total  ISR  for  each  case  is  30dB  with  a  processing 
gain  of  12dB.  The  signal  to  thermal  noise  ratio  is  varied  from  as  low  as  - 
lOdB  to  +5dB.  For  a  block  diagram  of  a  typical  system,  refer  to  [1]  with 
the  ARLO  detector  taking  the  place  of  the  LO  detector.  Also  included  are 
typical  spectra  for  the  signals  of  interest. 

In  general,  the  ARLO  detector  provides  robust  performance  in  all  in¬ 
terference  types.  It  outperforms  all  other  detectors  in  all  cases  with  the 
exception  of  the  single  CW  case,  where  it  closely  matches  the  performance 
of  the  LO  detector  without  memory.  In  the  partial  band  interference  sce¬ 
nario,  the  ARLO  detector  shows  drastic  improvement  over  the  LO  detector 
especially. 

7  Conclusions  and  Future  Work 

The  derivation  of  Eq.  (17)  provides  a  framework  for  simulation  of  the  autore¬ 
gressive  locally  optimum  nonlinearity  in  complex  interference  environments. 
Previous  simulation  [5]  focused  on  low  order  models  to  test  the  validity  of 
the  estimation  techniques.  As  can  be  seen  from  the  simulation  results  in 
Section  6,  the  performance  improvement  carries  over  into  more  complex  in¬ 
terference  scenarios.  As  anticipated  [3].  the  computational  complexity  of  the 
nonlinearity  is  significantly  reduced  from  an  iV -dimensional  jpdf  estimation 
to  a  P-dimensional  one. 

The  material  provided  in  this  final  report  represents  the  results  of  the 
work  performed  during  the  current  research  effort.  Given  these  results,  the 
following  issues  have  been  identified  as  possible  tasks  to  be  performed  in 
future  research  efforts: 
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Figure  1:  Partial  Band  environment  for  ISR  =  30dB  with  PG  =  12dB 
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Power  Spectra  for  PB  Interference  Environment 
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Probability  of  Bit  Error 


Detector  Performance  in  Single  CW  Interference  Environment 
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Figure  3:  CW  environment  for  the  case  of  ISR  =  30dB  with  PG 
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Figure  4:  CW  spectra  for  ISR  =  30dB  with  PG  —  12dB 
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Detector  Performance  in  Multi-CW  Interference  Environment 
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Figure  5:  Multi-CW  environment  for  total  ISR  =  30dB  with  PG  —  12dB 
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Power  Spectra  for  Multi-CW  Interference  Environment 
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Probability  of  Bit  Error 


Detector  Performance  in  Mixed  Interference  Environment 
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Figure  7:  Mixed  environment  for  ISR  =  30dB  with  PG  =  12dB 
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Power  Magnitude  (dB) 


Power  Spectra  for  Mixed  Interference  Environment 
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•  Task  1:  Simulation  of  other  interference  environments 

Although  several  scenarios  have  been  investigated  in  this  effort,  other 
interference  types  exist  that  require  study.  Examples  would  be  a  bursty 
jammer,  an  FM  modulated  jammer  and  groups  of  narrowband  jam¬ 
mers.  The  ARLO  detector  can  be  compared  against  the  LO  detection 
methods  of  previous  work  in  these  cases. 

•  Task  2:  Computational  Efficiency 

The  methods  of  computation  of  the  nonlinearity  have  been  improved 
by  the  translation  of  MATLAB  functions  to  C/C-b+.  To  implement 
LO  detection  within  an  embedded  system,  we  must  consider  the  issues 
of  memory  and  processing  speed.  Eventually,  the  algorithms  must  be 
performed  on  a  digital  signal  processor.  As  a  next  step  to  that  process, 
the  C  algorithm  development  should  be  continued  with  the  possibility 
of  implementation  on  a  DSP  and  real-time  performance  analysis  and 
feasibility  studies  should  be  performed. 

•  Task  3:  Application  of  other  LO  techniques 

The  ARLO  detector  has  been  developed  using  simple  histogram  pdf 
estimation  techniques.  Previous  work  [1]  suggests  other  means  of  esti¬ 
mating  the  nonlinearity  could  be  more  computationally  efficient  such  as 
the  FSA  method.  Enhancement  techniques  for  high  SNR  environments 
can  also  be  applied  to  the  ARLO  detector  to  increase  robustness.  For 
comparison,  performance  metrics  can  be  derived  for  the  ARLO  detector 
as  well. 

•  Task  4:  Projection  onto  convex  sets  (POCS) 

Projection  methods  have  been  used  to  detect  signals  in  narrowband 
interference  environments  with  promising  results.  In  particular,  the 
advantage  here  could  be  the  small  number  of  samples  needed  for  the 
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POCS  algorithm  in  estimating  the  pdf,  thus  reducing  the  real-time 
processing.  These  methods  should  be  investigated  in  conjunction  with 
the  ARLO  detection  scheme  in  other  interference  environments. 

•  Task  5:  Investigation  of  single  CW  interference 

The  ARLO  detector  does  not  perform  as  well  as  other  detectors  in 
the  single  CW  environment.  However,  it  outperforms  them  if  even  a 
second  CW  is  added  to  the  interference.  This  warrants  further  investi¬ 
gation  into  methods  to  improve  its  performance  in  the  single  CW  case. 
Possible  suggestions  include  varying  the  ”P”  parameter  in  coefficient 
estimation  techniques  or  changing  estimation  methods  altogether. 
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Appendix  1:  ARLO  Receiver  Block  Diagram 


Nonlinearity 


In  the  receiver  block  diagram  shown  above,  p  is  the  received  signal  vector 
of  length  N,  which  is  composed  of  both  the  signal  sent  by  the  transmitter 
and  the  noise  present  in  the  channel.  The  receiver  first  estimates  the  coeffi¬ 
cients,  {d},  in  the  auto-regressive  (AR)  model  using  the  modified  covariance 
algorithm.  Next,  the  AR  model  coefficients  and  the  received  signal  are  used 
to  estimate  h' ,  the  probability  density  function  (pdf)  of  the  noise.  This  is  cur¬ 
rently  accomplished  via  a  histogram  approximation  method.  The  resulting 
LO  nonlinearity,  g{p),  is  a  function  of  N  variables.  This  function  is  multiplied 
by  a  pseudonoise  (P-N)  sequence,  applied  to  a  summer,  and  then  a  decision 
is  made  as  to  whether  a  logical  ”0”  or  ”1”  was  sent. 
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Appendix  2:  MATLAB©  Simulation  Code 
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%  Program  silmulation .m 
%  August  15,  1998 
%  Myron  R.  Mychal 
% 

%  Simulates  AR  Model 

%  This  m-file  is  a  program  to  develop  a  robust  locally  optimum  receiver  for  a 
%  communications  system  modeled  in  a  colored  noise  channel.  This  program 
%  requires  the  functions  ar_histo.m  and  others.  The  non-linear 
%  processor  at  the  receiver  is  obtained  from  an  autoregressive  noise  model 
%  using  a  histcgram-based  estimation  technique  for  estimation  of  the  noise  pdf. 
% 

%  William  Jacklin,  Myron  Mychal 
%  Summer  '96  -  Summer  '97 

%  Set  parameters  for  simulation  run  -  initialization  procedure 

sim_parms;  %  Calls  up  the  required  parameters  to  run  the  simulation 

^  ★**★★*★*★*•*■★★★*★**★**★*★** 

%  **★★***★**★*****★*****★*** 

%  Begin  simulation  run  ** 

%  ************************** 

%  ****★★★***★*★*******★***** 

tic;  %  Start  stopwatch  to  time  one  run 

for  snr  =  1 : 1 : length (Eb_sigma)  %  Compute  results  for  all  Eb_sigma  values 
for  pe _ avg  =  1:1; it  %  Compute  average  of  (it)  iterations  of  Pe  for  each  snr 

%  ****************************** 

%  Implement  the  DSSS  transmitter 
%  ****************************** 


%  This  simulation  creates  a  BPSK  signal  of  "B"  bits  sampled  "Nb”  times 
%  for  each  iteration  of  the  BER  tester 

%  "ref"  is  the  BPSX  modulated  signal  of  randomly  created  information 

ref  =  2  *  rand ( 1 ,  B)  -  ones(l,  B) ;  %  Random  BPSK  data  for  1  iteration:  cts 

ref  =  sign(ref);  %  Data  symbols  at  1  sample  per  symbol 

sym_amp  =  sqrt(sigmat2  *  Eb^sigma { snr ) ) ;  %  Symbol  amplitude  assuming 

%  chip-matched  filtering 

%  IMPORTANT  NOTE:  Let  the  signal  have  amplitude  of  '1'  and  then  DIVIDE 
%  the  thermal  noise  by  the  symbol  amplitude  (below,  in  "Implement  the  Channel") 

[  coded^data,  chips  ]  =  dsss(ref,  Nc,  B,  PG,  1);  %  Transmitted  DSSS  signal 

%  and  PN  chip  sequence 


^  ********★★★★***★***★★ 
%  Ir.plement  the  channel 
^  *  *  ^*  ***************** 


%  Choose  the  interferer  type 
if  Itype  ==  'ar' 

[  interferer  ]  =  th_inter (Ntct ,  ISR,  a_t, 
elseif  Itype  ==  'pb' 

^  [  interferer  ]  =  pband(N,  ISR,  fc) ; 
elseif  Itype  ==  'ch' 

%  [  interferer  ]  =  chirp(N,  ISR,  fc) ; 
elseif  Itype  ==  'pi' 

%  ;  interferer  ]  =  pulsed (N,  ISR) ;  % 

elseif  Itype  ==  ' fx' 

%  [  interferer  ]  =  fmxnoise(N,  ISR);  % 
elseif  Itype  ==  'mx' 

[  interferer  ]  =  contwave (Ntot ,  ISRcw,  I, 
iidparms) ;  %  Mixed  interferers; 

elseif  Itype  ==  'cw' 


iidtype,  iidparms) ;  %  Theoretical  case 

%  Bandlimited  interference 
%  Chirped  interference 
Random  pulsed  interference 
FM  modulated  interference 

freq)  +  th_inter (Ntot ,  ISRpb,  a_t,  iidtype, 
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%  CW  interferers 


[  interferer  ]  =  contwave (Ntot ,  ISR,  I,  freq) ; 
end 

%  Add  thermal  noise  with  a  Gaussian  pdf  of  N(mu=0,  sigmat^2)  to  the  correlated 
%  interferer  (and  additive  channel  is  assumed) 

noise  =  (sigmat/sym_amp)  *randn (1 ,  Ntot)  +  interferer; 

%  implement  the  observed  signal 

%  ** -k -k  ****************** -k -kie  ** -k -k 

i  rcvd_tr  ]  =  receiver (coded_data,  noise,  T,  limit); 

^  ******************************************* 

%  Implement  a  linear  receiver  with  truncation 
^  *★**★★***★*****★*★★****★★★******★★★★★****** 

%  NO  AR  modeling  (whicener)  is  necessary 
%  g(rho)  =  rho,  so  no  histogram  is  necessary 

lc_dec  =  decision (chips ,  rcvd_tr,  Nb,  B) ;  %  Make  decisions 

num_errs_lc  =  sum  ( 0 . 5  *abs  ( lc_dec  ~  ref));  %  Find  no  of  errors  for  linear  receiver 
clear  lc_dec;  %  Memory  cleanup 

disp(['Done  with  linear  correlator.']); 

^  ******************************** 

%  Implement  LO  detector  techniques  %  THE  CRUX  of  this  whole  simulation 
^  ******************************** 

%  ?^############################ 

%  Find  the  aopropriate  whitener 

%  ##################ff#^######## 

%  We  should  loop  this  B/Bar  times 

a_hat  =  zerosd,  ?+l)  ;  %  Initialize  AR  parameter  estimate  vector 

[  a_hat  error  ]  =  modcov  (rcvd_tr  ( 1  :Ntot/ (B/Bar ))  ,  P)  ;  %  Find  the  whitener 

^  *★*******★*★★★*★**★★*★★★★******★*******★*★★★***★**★* 

%  Implement  a  linear  correlator  on  the  whitened  signal 
^  **************************************************** 

wc_dec  =  decision (chips ,  f ilter ( -a_hat ,  1,  rcvd_tr) ,  Nb,  B)  ;  %  Make  decisions 

num_errs_wc  =  sum(0 .  S’^abs  (wc_dec  -  ref));  %  Compute  number  of  errors 

clear  wc_dec;  %  Memory  cleanup 

disp(['Done  with  linear  correlator  on  whitened  signal']); 

g_lo  =  zerosd,  Ntot);  %  Initialize  the  LO  detector 

g_ar  =  g_lo;  %  Initialize  the  AR  detector 

g_lw  =  g_lo;  %  Initialize  the  LO  detector  for  the  filtered  case 

Nsamh  =  Bh  *  Nb;  %  Number  of  samples  in  each  histogram 

for  histonum  =  1:1:  (B/Bh)  %  Process  Bh*Nb  samples  at  a  time  for  a  histogram 

%  To  keep  track  of  which  iteration  the  simulation  is  on,  use  the  following  display 

%disp(['SNR  =  ',  num2str  (snr)  ,  '  It  =  '  ,  num2str  (pe^avg)  ,  '  Hist  =  ',  num2str  (histonum)  ]  ) 

%  ###################4#########^^ 

%  Implement  ARLO  and  LO  receiver 
%  ###################^^########## 

block  =  zeros (1,  Nsamh);  %  Initialize  block  to  prepare  for  the  histogram 

block  =  rcvd_tr(l,  ((histonum  -  l)*Nsamh+  1 ) : (histonum*Nsamh) ) ; %  Apply  the  correct  symbol 
[  h_lo,  bw_lo,  min_lo  ]  =  ar_histo (block,  Nb,  Bh,  -1,  K)  ;  %  LO  (P  =  0)  a_0  -  -1 

[  h_wh,  bw_wh,  rain_wh  ]  =  ar_his to  (block,  Nb,  Bh,  a_hat,  K)  ;  %  AR  methods 
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clear  block;  %  Memory  cleanup 

ivhitened  =  filter ( -a_hat ,  1,  rcvd_tr) ;  %  The  whitened  received  signal 
for  symcnt  =  l;Bh  %  Apply  nonlinearity  to  each  symbol 

symbol  =  zeros (1,  Nb) ;  %  Look  at  one  symbol 

symbolindx  =  ( (histonum-l) *Nsamh+ { symcnr-1) *Nb) +1  :  ( (his tonum-1 ) *Nsamh+ ( symcnt ) *Nb) 

symoolindx  =  symbolindx ;  %  Force  it  to  be  a  row  vector 

symbol  =  rcvd_tr ( symbolindx) ;  %  Get  the  current  symbol 

symbol  =  symbol ;  %  Force  it  to  be  a  row  vector 

filsym  =  whitened ( symbolindx) ;  %  The  whitened  symbol 

filsym  =  filsym(:)';  %  Force  it  to  be  a  row  vector 

g_lo(l,  symbolindx)  =  nlp(symbol,  min__lo,  Nb,  h_lo,  bw_lo,  -*1)  ; 

g_lw(l,  symbolindx)  =  nlp(filsym,  min_wh,  Nb,  h_wh,  bw_wh,  -1) ; 

g_ar{l,  symbolindx)  =  nip (symbol,  min_wh,  Nb,  h_wh,  bw_wh,  a_hat) ; 

end 

end  %  End  of  loop  on  "numhistos ” 

%  ##^################################ 

%  Make  ARLO  and  LO  detector  decisions 
%  ##??##############  =  ################# 

lo_dec  =  decision (chips,  g_lo,  Nb,  B) ;  %  Make  decisions  for  LO 

num_errs_lo  =  sum( 0 . 5*abs (lo_dec  -  ref));  %  Count  no  of  errors  for  LO 

clear  g_lo,  lo_dec;  %  Memory  cleanup 

lw_dec  =  decision (chips,  g__lw,  Nb,  B) ;  %  Make  decisions  for  LO  on  whitened 

num_errs_lw  =  sum(0 . 5*abs (lw_dec  -  ref));  %  Count  no  of  errors  for  LO  on  whitened 

clear  g_lw,  lw_dec;  %  Memoiry  cleanup 

ar_dec  =  decision (chips ,  g_ar,  Nb,  B) ;  %  Make  decisions  for  ARLO 

num_errs_ar  =  sum. ( 0 . 5*abs (ar_dec  -  ref));  %  Count  no  of  errors  for  ARLO 

clear  g_ar,  ar_dec;  %  Memory  cleanup 

%  Compute  probability  of  bit  error  for  this  iteration  and  this  "snr" 

Pb_ic  (pe_avg,  snr)  =  num__errs__lc  /  B;  %  Linear  correlator 

?b_v/c  (pe_avg,  snr)  =  num_errs_wc  /  B;  %  Linear  correlator  on  whitened 

Pb_lo (pe_avg,  snr)  =  num_errs__lo  /  B;  %  LO  memoryless  detector 

?b_lw ( pe__avg ,  snr)  =  num_errs_lw  /  B;  %  LO  detector  on  whitened  signal 

Pb_ar (pe^avg,  snr)  =  num_errs_ar  /  B;  %  ARLO  detector  with  memory 

end  %  End  of  iteration  on  "pe_avg” 

end  %  End  of  iteration  on  ’‘snr" 


%  Compute  average  Pb  (from  Monte  Carlo  analysis) 


if  (it  ==  1) 
?b_lc_avg  = 
Pb_v/c_avg  = 
Pb_lo_avg  = 
Pb_lw_avg  = 
Pb_ar_avg  = 
else 

Pb_lc_avg  = 
Pb_wc_avg  = 
?b_lo_avg  = 
?b_lw_avg  = 
?D_ar_avg  = 
end 


%  If  only  one  iteration 

Pb_lc ; 

Pb_wc ; 

Pb_lo; 

Pb_lw ; 

Pb_ar ; 

%  If  several  iterations 
(1/it) *sum ( Pb_lc) ; 

(1/it)  ’^sum.(Pb_wc)  ; 

(1/it)  ’‘^sum( Pb_lo)  ; 

(1/it) *sum(Pb_lw) ; 

(1/it) *sum(Pb_ar) ; 


run_cime  =  toe;  %  Stop  stopwatch 

save  cwr08  Pb_lc_avg  Pb_wc_avg  Pb_lo_avg  Pb_lw_avg  Pb_ar_avg  a_hat  run_time; 
%save  ewsOS  noise  rcvd_tr  interferer; 


%  PLOT  resulting  Pb  curves 


%f igure; 
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%hold  on; 

%semilogy ( Pb_lc_avg, 

%tit:le  (' Detector  Performance  in  INSERT  Interference  Environment'); 

%xlabel ( ' Signal-to-Thermal  Noise  Ratio  (dB)'); 

%ylabel (' Probability  of  Bit  Error'); 

% semi logy (Pb_wc_avg,  'c-'); 

%semilogy (Pb_lo_avg,  'r-' ) ; 

%semilogy (Pb_lw_avg,  'w-'); 

%semilogy (Pb_ar_avg,  ' g- ' ) ; 

%legend{gca,  'Linear  Correlator',  'Correlator  on  whitened',  ' LO  detector ' ,  'LO  on 
whitened' ,  'ARLO' ' ; 

%hold  off; 

end 
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%  This  is  a  data  file  that  contains  the  simulation  parameters  for  the  ARLO  detector 
%  The  file  "simulation. m"  calls  this  file  and  some  of  the  functions  within 
%  " simulation .m"  require  these  global  definitions.  See  the  individual  functions 
%  for  comments  on  where  each  parameter  is  used 
%  Book-keeping  parameters 

fnair.e  =  'ARresults';  %  Output  filename 

it  =  1;  %  Number  of  Monte  Carlo  iterations  to  average  for  Pe  calculation 

%  pdf  estimation  technique  parameters 
%  -  also  see  ar_histo.m,  lo_histo.m  &  receiver. m 

Bh  =  512;  %  Total  number  of  symbols  to  calculate  histogram 

K  =  33;  %  Data  points  to  calulate  (#  of  breakpoints) 

T  =  1000;  %  Range  of  support  for  pdf 

%  Spread-spectrtim  parameters 

Nc  =  4;  %  Samples  per  chip 

PG  =16;  %  Processing  gain  (chips  per  bit) 

Nb  =  PG*Nc;  %  Samples  per  symbol  (assume  symbol  is  a  bit) 

B  =  Bh;  %  Symbols  transmitted/received  and  sent  to  BER  tester 

Ntot  =  B  *  Nb;  %  Total  number  of  samples  in  run 

%  AR  model  parameters 

Ear  =  Bh/64;  %  Total  number  of  symbols  to  estimate  AR  parameters 

Bg  =  Bh;  %  Total  number  of  symbols  to  calculate  nip 

ar_method  =  'mcov' ;  %  Coefficient  estimation  method 

%  -  'mcov'  =  Modified  Covariance 
%  -  'cov'  =  Covariance 
%  -  'yw'  =  Yule  Walker 

%  -  'burg'  =  Burg  Method 

%  Signal  parameters  and  thermal  noise 
%  -  also  see  receiver. m 

sigmat  =  1;  %  Standard  deviation  of  thermal  noise 

sigmat2  -  sigmat''2;  %  Variance  of  thermal  noise 

Eb_min  =  -10;  %  Minimum  value  of  thermal  SNR  in  db 

Eb__max  =0;  %  Maximum  value  of  thermal  SNR  in  db 

Eb_sigma  =  Eb_min : Eb_max ; 

for  index  =  1 : length (Eb_sigma) 

Eb_sigma( index)  =  10-"  (Eb^sigma  (index) /lO)  ;  %  Vector  of  SNR's 

end 

limit  =  'hard';  %  Limiter  on  noise 

%  -  'hard'  =  noise  blanker 
%  -  'soft'  =  soft  limiter 

clear  index; 

%  Correlated  interferer  parameters 

%  Interf erer-to-signal  ratio  (db) 

%  Interference  type 
'cw'  =  Continuous  Wave 
'ar'  =  Theoretical  AR 
'pi'  =  Pulsed 
'pb'  =  Partial  Band 
'ch'  =  Chirped  FM 
'fx'  =  FM  Modulated  Noise 

%  'mx'  =  mixed  partial  band  and  CW 

%  Continuous  wave  jammers 

'mx'  ) 

%  Number  of  CW  interferers 
%  CW  Interferer  frequencies 
%  Number  of  AR  coefficients 
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if  (Itype  ==  'cw')  |  (Itype  == 

1  =  1; 

freq  =  [  Nc*B  ]  ; 

P  =  16; 


ISR  =  30; 
Itype  =  'cw' 

% 

% 

% 

% 

% 

% 


ena 


%  Wide-band  filtered  white-noise  jammers 
if  (Itype  ==  'ar')  |  (Itype  ==  'mx') 

poly  =  [  1  -0.75];  %  The  pole  of  the  low-pass  filter 

L  =  4;  %  For  a  filter  of  order  2^L 

for  i  =  1:L 

poly  =  conv(poly,  poly);  %  Recursively  increase  the  order  of  the  filter  to  2^L 
end 

clear  i; 
a_t  =  -poly; 

P  =  length(a_t)  -  1;  %  The  number  of  coefficients 

iidtype  =  'gs';  %  The  type  of  iid  noise  sequence 

%  -  'gs'  =  Gaussian 

%  -  'pp*  =  Poisson 

%  -  'Ip'  =  Laplacian 

%  -  'un'  =  Uniform 

%  -  'mg'  =  Mixed  Gaussian 

%  -  'ey'  =  Cauchy 

%  -  'rl'  =  Rayleigh 

%  -  'ri'  =  Ricean 

%  -  'nk'  =  Nakagami-m 

iidparms  =[01];  %  Parameter  vector  for  iidtype 

%  -  for  'gs'  =  [  mean  standarddev  ] 

%  -  for  'pp'  =  [  rate_parameter  ] 

%  -  for  'Ip'  =  [  sigma  ] 

%  -  for  'un'  =  [  mean  length  ] 

%  -  for  'mg'  =  [  meanl  stddevl  mean2  stddev2  ] 

%  -  for  'ey'  =  [  mean,  scale_parameter  ] 

%  -  for  'rl'  =  [  ] 

%  -  for  'ri'  =  [  ] 

%  -  for  'nk'  =  [  ] 

end 

%  Mixed  jammers 
if  (Itype  ==  'mx' ) 

ISRcw  =  ISR/2; 

ISRpb  =  ISR/2; 

end 
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%  This  M-file  develops  the  decision  statistic  for  the  various  types  of 
%  receivers  in  a  communication  system.  It  assumes  that  some  function 
%  of  the  received  signal,  g{rho),  is  used  in  a  direct  sequence  spread 
%  spectrum  environment  and  multiplied  by  a  chip  sequence.  This  is  then 
%  passed  through  a  discrete  integrator  (a  summer)  and  sampled  at  the 
%  conclusion  of  the  symbol  (after  Nb  samples  have  passed) .  The  decision 
%  assumes  synchronous  BPSK  modulation. 

% 

%  Inputs : 

% 

%  pn  -  The  PN  chip  sequence 

%  g  -  The  detector  function,  g(rho) 

%  Nb  -*  The  number  of  sample  points  per  symbol 

%  B  -  The  number  of  symbols 

% 

%  Outputs : 

% 

%  dec  -  The  decision  vector 

function  [  dec  ]  =  decision (pn,  g,  Nb,  B) 

dec  =  ones(B  ,1);  %  Assume  that  all  symbols  are  "I's” 

stat  =  sum ( reshape  (pn  .*  g,  Nb,  B)  )  '  ;  %  Sum  each  symbol  correlated  with  the  pn 

HO  =  find(stat  <=  0);  %  ML  estimation  choosing  where  "-I's"  are 

dec(HO)  =  -ones (size (HO) ) ;  %  CHange  symbols  accordingly 

dec  =  dec(:)';  %  Make  it  a  row  vector 

end; 
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Appendix  3:  C  Simulation  Code 
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Appendix 


/*  filter. c 

*  Version  1.0,  3/22/98 

*  Written  by  Michael  Banys  &  Angela  Kaczmarski 

* 

*  Revision  Control: 

* 

* 

* 

*/ 


/★  ★★******★*»★*★***★***★**★★★★******★***★**★★*★★**★★★**★***★/ 

*  FILTER  Function 

*  OUTPUT  =  filter {  b[M] ,  a[N].  data_in[]  ) 

* 

*  b[M]  =  coefficients  of  numerator  of  impulse  response 

*  atN]  =  coefficients  of  denominator  of  impulse  response 

*  =  order  of  filter 

* 

*  This  function  computes  the  convolution 

*  of  vector  '’data_in"  with  the  causal  filter 

*  whose  time-domain  filter  coefficients 

*  are  of  the  following  format: 

* 

*  a(l)*y(n)  =  -SUM(k=l,N)  [  a(k)*y(n-k)  ] 

*  +SUM(k=0,M)  [  b(k)*x(n-k)  ] 

* 

*  If  a{0)  !=  1,  FILTER  will  normalize  the  filter 

*  coefficients  by  a(0) .  If  a(0)==0,  return  ERROR. 

•k 

*/ 


#define  NULL  'OL' 

#define  ORDER_ERROR  "Filter  numerator  order  is  greater  than  denominator 

order\n" 

double  filter (  double  b[],  double  a[],  double  data_in[],  double  y[]) 

{ 


int  f ilter_order , 

f ilter_iteration=0 , 
element_num, 


i  =  0, 

M, 

N, 

INPUT_LENGTH; 


/*  Order  of  acting  filter  */ 

/*  Variable  tracking  loop  through  input  data  */ 
/*  Variable  indicating  current  coefficient  and 
delay  being  processed  */ 

/*  Iteration  variable  */ 

/*  Order  of  Numerator  */ 

/*  Order  of  Filter  (Denominator)  */ 

/*  Length  of  input  data  vector  */ 


double  normalization^value,  /*  a[0]  coefficient  v^lue  */ 
*f ilter_delay,  /*  Filter  delay  values  */ 
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sumx, 


/*  Temporary  variable  representing  sum  of 
elements  of  a[]  multiplied  by  elements  of 
f ilter_delay [ ]  */ 

sumy;  /*  Temporary  variable  reperesenting  sum 

of  elements  of  b[]  multiplied  by  elements 
of  f ilter_delay [ ]  */ 


/*  Determine  number  of  elements  in  array  b[]  */ 
for  (i=0;  b[i]  !=  NULL;  i++) 

t 

M  =  i; 

/*  Determine  number  of  elements  in  array  a[]  */ 
for  (i=0;  a[i]  !=  NULL;  i++) 

I 

N  =  i; 

/*  Determine  number  of  elements  in  array  data_in[]  */ 
for  (i=0;  data_in[i]  !=  NULL;  i++) 

/ 

INPUT_LENGTH  =  i; 

filter_delay  =  malloc (N  *  sizeof (double)  );  /*  Allocate  memory  for 

filter_delay  length  */ 

filter_order  =  N;  /*  Set  order  of  filter  */ 

/*  Check  to  see  if  the  length  of  b[]  <  length  of  a[] 

If  length  b[]  >=  length  a[],  output  error  message  and  exit  function  */ 
if  (M  >=  N) 

{ 

printf (ORDER^ERROR) ; 
return  0; 


/*  Check  to  see  if  denominator  of  IR  is  normalized  (i.e.  does  a[0]=l 
If  denominator  is  not  normalized,  normalize  IR  by  a[0]  */ 
if  (a[0] !=  1) 

{  normalization_value=a [0] ; 

/*  Normalize  denominator  */ 
for  (i=0;  i<N;  i++) 

a [i] =a [i] / norma lization_value; 

/*  Normalize  numerator  */ 
for  (i=0;  i<M;  i++) 

b [ i ] =b [ i] /normalization_value; 


/*  Set  w[n]  delay  array,  f ilter_delay [ ] ,  to  all  zeroes  */ 
for  (i=0;  i<=N;  i++) 
f ilter_delay [i] =0 ; 
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/*  start  filtering  process  */ 

for  ( f ilter_iteration=0;  filter_iteracion<INPUT_LENGTH; 
f ilcer_iteration++) 

{ 

sumx=sumy=0.0;  /*  Initialze  variables  */ 

/*  Sum  left  side  of  Direct  Form  II  Structure  */ 
for  ( element_num=:l ;  element_num<N;  element_num+  +  ) 

sumx  =  sumx  -  a  [element^num]  *  f  ilter^de lay  [element _num]  ; 

/*  Sum  right  side  of  Direct  Form  II  Structure  */ 
for {element_num=l;  elemenc_num<M;  element_num++) 

sumy  =  sumy  +  b  [element_num]  *  f ilter_delay  [element_num]  ; 

/*  Compute  current  filter  delay  {f ilter^delay [0] )  */ 
filter_delay [0]  =  data_in [ f ilter_iteration]  +  sumx; 

/*  Compute  output  element,  y(n)  (y [ filter_iteration] )  */ 
y[filter_iterationl  =  b[0]  *  f ilter^delay [ 0]  +  sumy; 

/*  Delay  each  f ilter_delay [ ]  element  by  one  */ 

for  (  element„num=N-l;  element_num>0 ;  element__num —  ) 

f  ilter_delay  [element_num]  =f  ilter__delay  [element_num-l ]  ; 


} 

free ( f ilter_delay)  ; 
return  0; 


} 
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/*  hist.c 

*  Version  1.0,  8/24/98 

*  Written  by  Michael  Banys  &  Angela  Kaczmarski 

* 

*  Revision  Control: 

* 

* 

* 

*  / 

/ 


/*  hist(  double  signal[],  int  bin^info,  double  f req_count [ ] ,  int  bin_loc[]  ) 
* 


*  freq__counr:  =  vector  containing  frequency  distribution  of  each  bin 

*  bin_loc  =  vector  containing  center  location  of  each  bin 

* 

*  signal  =  signal  (vector)  from  which  to  compute  histogram 

*  bin_info  =  number  of  bins  necessary  for  histogram 

* 

*  NOTE:  All  vectors  must  be  NULL  terminated  for  size  determination 

★ 

*  The  histogram  function  bins  the  elements  of  signal  into  the  number  of  bins 

*  specified  by  bin_info.  The  function  returns  the  number  of  elements  in 

*  each  bin  and  the  center  point  of  each  bin. 

*/ 


#include  <stdio.h> 
#include  <stdlib.h> 

#define  NULL  'OL' 


void  hist(  double  signal [], 
{ 

int  signal_length, 

i, 

num_bins , 

j  / 

accumulator, 

* temp_f req_count ; 


int  bin_info,  int  f  req_^count  [  ]  ,  double  bin_loc  [  ]  ) 

/*  length  of  signal  vector  */ 

/*  counting  variable  */ 

/*  Number  of  bins  */ 

/*  Intermediate  bin  center  locations  */ 


double  min_y, 
max_y, 

temp_double, 
*temp_bin_loc , 
binwidth; 


/*  Minimum  signal  value  */ 
/*  Maximum  signal  value  */ 

/*  Width  of  each  bin  */ 


for  {i=0;  signal(i]  !=  NULL;  i++) ;  /*  Determine  length  of  signal[]  */ 

signal_length  =  i; 

min_y  =  signalCO];  /*  Initialize  min_y  */ 

for  (i=l;  i  <  signal_length;  i++)  /*  Find  actual  min_y  */ 

( 

temp_double  =  signal [i ] ; 
if  (temp_double  <  min  v) 
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} 


inin_y  =  temp_double; 


max_y  =  signal[0];  /*  Initialize  inax_j/  */ 

for  (i=l;  i  <  signal_length;  i++  )  /*  Find  actual  may:_y  */ 

{ 

teinp_double  =  signal  [ i]; 
if  ( temp_double  >  max^) 
max  _y  =  temp_double; 

} 

binwidth  =  (max_y  -  min_y)  /  bin_info;  /*  Determine  binwidth  */ 

temp_bin_loc  =  malloc (  (bin_info+l)  *  sizeof (double)  ); 

/*  Allocate  memory  for 
temp_bin_loc  */ 

for  {i=0;  i  <=  bin_info;  i++)  /*  Determine  bin  centers  */ 

temp_bin_loc[i]  =  min_:y  +  binwidth*i; 
temp_bin_loc [ (int) bin_info]  =  max_y; 

for  (i=0;  i  <  bin_info;  i++) 

bin_loc[i]  =  temp_bin_loc [i]  +  (binwidth/  2.0); 


temp_freq_count  =  malloc(  (bin_info+l)  *  sizeof (int)); 


for(  i=0;  i  <  signal_length;  i++) 
temp_frecucount [i]  =  0; 


/*  Fill  first  row  of 
temp_freg_count  with  zeros  */ 


for(  i=l;  i<bin_info+l;  i++)  /*  Determine  bin  member  info  */ 

{ 

accumulator  =0;  /*  Initialize  accumulator  */ 

for(  j=0;  j  <  signal_length;  j++) 

if(  signal [j]  <=  temp_bin_loc [i]  ) 

accumulator++ ; 

temp_freq_count [i]  =  accumulator; 

} 

for  (i=0;  i<bin_info;  i++)  /*  Determine  exact  bin  member  info  */ 

f req^counc [ i ] =temp_f req_count [ i+1 ] - temp_f req_count [ i ] ; 


free ( temp_freq_count) ;  /*  Free  dynamically  allocated  memory  */ 

free(temp_bin_loc) ; 


} 
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Appendix  4:  C++  Simulation  Code 
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'  Crsaced  By  Abhijic  Talachi  for  £C5  449 
Simuiation  of  DSS3  Communica cions  Syscem 

This  header  file  contains  the  class  definition  for  an  Interferer 
The  interferer  is  derived  from  the  class  Signal 
Three  types  of  interferers  are  available. 

1.  Gaussian 

2 .  Uniform 

2 .  Continuous  Wave 

iifndef  INTERFERER_H 
^define  INTERFERER^H 

4  include" signal .h" 

const  int  MAX_FREQ  =  64000; 
const  int  MAX_LEVEL  =  100; 

cempiate<class  T>  class  Interferer  :  public  Signal<T> 
protected: 

//  Interference  to  signal  ratio 

double  isr;  //  Interference  to  signal  ratio 

//static  const  int  MAX_FREQ  =  64000;  //this  is  not  accepted  by  complier 
//static  const  int  MAX_LEVEL  =  100;  //so  is  globally  defined 

public ; 

//  Constructor 

Interferer { cons t  long  length,  const  double  level  =  1); 


;/  Constructor  for  interferer 
template<class  T>  Interf erer<T> 


Interferer {const  long  length,  const  double  level ); Signal<T> (length) 


try 

{ 

if (level  >  MAX_LEVEL  ) 
throw  " error " ; 

} 

catch  (...) 

^  cout«endl« “Interferer  :  Bad  argximent  for  interferer  (ISR  is  too  large*  )“; 
exit  (EXIT_F.\ILURE)  ; 


level ; 


'/  Continuous  wave  Interferer 

cemplate<class  T>  class  Contwave  :  public  Interf erer<T> 

{ 

public : 

//  Constructor  ^  ^  r  \ 

Contwave { cons t  long  legth,  const  double  isr,  const  douDle  freg) , 


/ /Constructor 


//  Function  does  nothing  for  CW  interferer  as  attenuation  is  accoxinted  for 

//  in  the  constructor 

void  attenuate (const  double  factor) {} 


//  Overloaded  >>  operator  for  user  derined  CW  interferer 
friend  istream  &operator  >>  (istream  ies tream, Contwave<T>  Scproto)  ; 


■ /  Constructor  for  the  class  Interferer 

•  /  This  irav  be  modified  later  to  take  a  number  of  frequencies  and  phases  ^  ^  ,  -r 

tempiate<class  T>  Contwave<T> :: Contwave (const  long  length,  const  double  leve^,  const  double  freq):^n 

terf erer<T> ( length, level ) 

{ 

double  phase  =  0.0,  time,  amp; 

//  Phase  will  be  made  a  vector  of  phases  if  required 

amo  =  sqrt(2  *  pow(10,  (isr  /  10))); 

double  deita^t  =1.0/  length;  ;!  Time  interval 

for  (int  i  =  0;  i  <  length;  1-^+) 


daca[i]  =  amp  *  sin  (2  *  PI  *  freq  "  time  +  phase);//  Sinusoidal  interferer 
Interferer" ) ; 


name("Cont  v;ave 
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Q  n 


Contwave<T>  &o) 


Overloaded  >>  operacor  for  user  defined  CW  incerferer 
“emplace  <class  T>  iscream  Scoperaccr  >>  (iscream  Scstream, 


long  lengthi- 
double  level, freq; 


couC<< “Creating  a  Continuous  wave  interferer*  <endl; 

cout<<‘' (This  may  be  modified  later  for  multi':. e  interferers)  “ <<endl ; 

cout<<“ Enter  the  Number  of  samples “<<endl; 
s  t r eam>  > 1 ength ; 

cout<<“Encer  the  Interference  to  Signal  Ratio  (ISR) ”<<endl; 
stream>> level ; 


cout<< “Enter  the  frequency  of  the  wave“<<endl; 
stream>>f req; 


try 


if  (level  >  MAX_LEVEL  |  |  freq  >  MAX^FREQ  ) 
throw  “  error  *•  ; 

/ 

catch ( . . . ) 

{ 

couc<<endl«"ContWave  Inter:  Bad  arguments  for  CW  interferer"; 
exit (EXIT^FAILURE) ; 


Contwave<T>  Temp (length, level , freq) ; 

//Fill  in  the  elements  of  the  continuous  wave  interferer. 
o=:Temp; 


return  stream; 

} 

’/  Theoretical  Gaussian  Interferer 

tempiate<class  T>  class  Gaussinter  :  public  Interf erer<T> 

{ 

public: 

Gaussinter (const  long  legth,  const  double  isr,  const  double  mu  =  0.0,  \ 
const  double  sigma  =1.0); 

friend  istream  &operator  >>  (istream  Scstream,  Gaussinter<T>  Stproto)  ; 

//  Overloaded  >>  operator  for  user  defined  Gaussian  interferer 

\  ; 

emplate<class  T>  Gauss inter<T> :: \ 

aussinter (const  long  length,  const  double  level,  const  double  mu  ,  const  double  sigma) \ 
:  Interf erer<T> ( length, level ) 

gaussgen ( length,  mu,  sigma  ) ; 

attenuate (isr) ; 

name ( "Gaussian  Interferer" ) ; 


template<ciass  T>istream  &operator>> ( istream  ^stream,  Gaussinter<T>  &o) 
long  le.ngth; 

double  level,  mu,  sigma; 

cout<< “Creating  a  Gaussian  Noise  Sequence " <<endl ; 

cout<< “Enter  the  Number  samples “ <<endl ; 

3  tr eam>> 1 ength ; 

cout<<"Enter  the  Interference  to  Signal  Ratio  ( ISR) " <<endl ; 
stream>> level ; 

couC<<“ Enter  the  Mean  of  the  linterf erer " <<endi ; 
s  tream>>mu; 


cout<<" Enter  the  Standard  Deviation  of  the  Interf erer "<<endl ; 
stream>>sigma; 

..  /  Exceptions  for  arguments  are  handled  in  the  function  GaussgenO 
Gaussincer<T>  Temp ( length, level , mu , sigma) ; 

'/Fill  in  the  elements  of  the  Gaussian. 
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o=Temp; 

return  stream; 

} 

/  Theoretical  uniform  interference 
tempiate<ciass  T>  class  Unifointer  :  public  Interf erer<T> 

{ 

Unifointer (const  long  legth,  const  double  isr,  const  double  mu  =  0.0); 

//  Overloaded  >>  operator  for  user  defined  Uniform  Interf erer 
friend  istream  Scoperator  >>  (istream  Scstream,  Unifointer<T>  &proto)  ; 


} ; 


template<class  T>  Unifointer<T> 


rUnifointer (const  long  length,  const  double  level,  \ 

const  double  mu  ): Interf erer<T> ( length, level) 


{ 

unifcgen ( length,  mu) ; 

attenuate ( isr)  ; 

name ( “Uni form  Interf erer “ ) ; 

} 


//  Overloaded  istream  operator  for  an  object  or  typ©  Unifointer 
template<class  T>  istream  &operator»  ( istream  ScStream,  Unif ointer<T>  &o) 

{ 

long  length; 
double  mean, level; 


cout<<  “Creating  a  Uniform  Interf  erer '‘<<endl  ; 

cout«“ Enter  the  Number  of  samples  "«endl  ; 
s  tream>> 1 ength ; 

cout«“ Enter  the  ISR"«endl; 
stream»level  ; 

cout«"  Enter  the  Mean  the  Interf  erer  “«endl  ; 
stream»mean; 

//  Exceptions  for  argximents  are  handled  in  the  function  UniforgenO 
/ /  Create  an  object  of  type  Uniform  Interference 

//  Exceptions  for  arguments  are  handled  in  the  function  Gaussgen() 
Unif ointer<T>  Temp (length, level , mean) ; 

//Fill  in  the  elements  of  the  continuous  wave  interf  erer . 
o=Temp ; 

return  stream; 


#endif 
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treated  By  Abhijit  Taiathi  :cr  BCE  449 
Simula:: ion  of  DSS3  Communicacicns  System 

This  header  file  contains  the  class  definition  for  Noise 
The  Noise  is  derived  from  the  class  Signal 
Tv'pes  of  Noise, 

1.  Gaussian 

2.  Uniform 


Tifndef  !I0I3E_H 
4 define  NOISE^H 


4 include "signal .h" 


Class  for  Noise 

tempiate<ciass  T>  class  Noise  :  public  Signal<T> 
protected: 

long  seed;  //  Seed  for  Noise  generation 
public : 

Noise(const  long  length); 


template<class  T>  Noise<T> : :Noise (const  long  length) : Signal<T> (length) 

{ 

seed  =:  (int)  time  (NULL)  ; 
srand(seed)  ; 

} 

//^Gaussian  Noise 

cemDiate<class  T>  class  Gaussnoise  :  ciiblic  Noise<T> 

{ 

public : 

Gaussnoise (const  long  legth,  const  double  mu  =  0.0,  const  double  sigma  =  1.0); 

friend  istream  &operator  >>  (istream  &stream,  Gaussnoise<T>  &proto) ; 
//Overloaded  >>  operator  for  user  defined  Gaussian  Noise 


template<class  T>  Gaussnoise<T> :: Gaussnoise (const  long  length,  const  double  mu,  const  double  sigma) 
Noise<T> (length) 

gaussgen( length,  mu,  sigma  ); 
name ( “Gaussian  Noise" ) ; 

template  <class  T>istream  &operator>>  ( istream  istream,  Gaussnoise<T>  &o) 

{ 

long  length; 

double  level,  mu,  sigma; 

cout« “Creating  a  Gaussian  Noise  Sequence "<<endl; 

cout<<‘ Enter  the  Number  samples " <<endl ; 
stream>>length; 

ccut<< “Enter  the  Mean  the  Noise" <<andl ; 

3tream>>mu ; 

cout<<  ■' Enter  the  Standard  Deviation  of  the  Noise  *' <<endl  ; 

3 tream>>sigma ; 

Exceptions  are  handled  in  the  gaussgen  function 
Gaussnoise<T>  Temp  ( length, .mu,  sigma)  ; 

/■Fill  in  the  elements  of  the  continuous  wave  interferer. 
c  =Temp ; 

return  stream; 


'.Uniform  noise 

tempiate<ciass  r>  class  Unifonoise  :  public  Noise<T> 

( 

cubl ic : 

Unifonoise (const  long  iegth,  const  double  mu  =  0.0); 

friend  istream  ioperator  >>  (istream  istream,  Unifonoise<T>  iproto); 
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//Overicaded  >>  operacor  for  user  defined  Uniform  Noise 


;emplate<class  T>  Unif onoise<T> : : unif onoise (consr  long  length,  const  double  mu) :Noise<T> (length) 

unifogen { length,  mu); 
name { -Uniform  Noise”); 

} 

//Overloaded  iscream  operator  for  an  object  of  type  Unifonoise 
template<class  T>istream  Scoperator»  ( istream  Scstream,  Unif onoise<T>  &o) 

{ 

long  length; 
double  mean; 

oout« "Creating  a  Uniform  Noise  Sequence “<<endl; 

cout<<  **  Enter  the  Number  of  samples ''<<endl ; 
stream>> length; 

cout«'' Enter  the  Mean  the  Noise” <<endl ; 

3tream»mean; 

Unifonoise<T>  Temp (length, freq) ;  ^ 

//Fill  in  the  elements  of  the  continuous  wave  interferer. 

o=Temp; 

/ /  Exceptions  are  handled  in  the  Gaussgen  runction 
return  stream; 

} 

?^endif 
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/ 

■  / 


Created  3y  Abhijit  Talathi  for  £CE  449 
Simulation  of  D3S3  Communications  System 

This  header  file  contains  the  class  definition  for  a  BPSK  signal 
The  BPSK  is  derived  from  the  class  Signal 
This  cosists  of  a  binary  signal  composed  of  the  symbols  -1  and  1 

'/  Class  binar-'/  chase  shift  keyed  signal  {BPSK)  (the  message  encoded  in  ones  and  -ones  ) 
#ifndef  3PSK_H 
define  BPSK_H 

# include "signal .h“ 

template<class  T>  class  BPSK  :  public  Signal<T> 
public : 

//  Constructor 
BPSK(const  long  length); 

//  Generates  a  BPSK  signal 
void  BPSKgen (const  long  length) ; 

//  Performs  Direct  Sequence  Spread  Spectrum  codification  of  the  BPSK  signal 
Signal<T>  dsss(  Signal<T>  ^signal,  const  int  KC, 

const  int  B,  const  int  PG  ,  const  int  amp) ; 

//  Overloaded  »  operator  for  user  defined  BPSK  signal 
friend  istream  ^operator  >>  {istream  Scstream,  BPSK  &proto) ; 


template<class  T>  3PSK<T>  ::  B?SK(const  long  length) : 3ignal<T> ( length) 
{ 

BPSKgen (length) ; 
name ( "BPSK* ) ; 


tsmpiate<class  T>  void  3PSK<T>  ::  BPSKgen ( const  long  length) 

for (int  i=0;  i  <  length;  i++) 

{ 

data[i]  =  (rand()%2)  ?  1  :  -1;//  Initailize  to  ones  and  -ones 

} 


Implement  the  function  DSSS  ie  Discrete  Spread  Spectrum  Signal  on  The  BPSK  signal 
template<class  T> 

5ignai<T>  3PSK<T> : : dsss (  Signal<T>  isignal . const  int  NC, 

const  int  3,  const  int  PG  ,  const  int  amp) 

ConvertToRowMatrix ( ) ;  //  Method  of  class  Matrix 

int  MB  =  NC  PG;  ,  i  Constants  for  DSSS 
int  N  =  NB  »  3; 

Matrix<T>  tmp (NB , 1 , 1 ) , mat 1 ; 

tmp=tmp  *(*this);  Multiply  the  Matrices 

matl=tmp; 

matl .ConvertToRowMatrix ( ) ; 

tmp  =  RandomMatr ixOne  ( 1 ,  ( 3*  PG )  }  "^  2  ; 

tmp-=l;  //  Subsract  1  from  temp  to  get  BPSK 

tmp.SignO;  //  Final  BPSK  signal 

tmp=  Matrix<T> (NC, 1 , 1} *tmp;  /  '  Multiply  by  ones 

Matrix<T>  chips=tmp; 

chips . Reshape iN. 1 ) ;  /,  Invoke  reshape  in  the  parent 
chips . Sign () ; 

chips . ConvertToRowMatrix ( ) ; 

Matrix<T>  Temp  =  chips;  Create  signal 
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Temp . Indmui t ( mac 1 ) ; 


3ignai=Temp  *  amp; 
return  chips; 

} 

//  Overloaded  iscream  operator  for  an  ooject  of  type  BPSK 
template<class  T>iscream  &operator» ( is cream  ^stream,  5PSK<T>  &o) 
{ 

long  length; 

cout« "Creating  a  Binary  Phase  Shift  Keyed  Sequence " «endl ; 

cout« "Enter  the  Number  of  samples‘'<<endl; 

3  cream»  length; 

//  Create  an  object  of  type  BPSK 
o.BPSKgen( length) ; 

return  stream; 

} 

#endif 
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constructs  the  final  received  signal 
\piate<ciass  •?> 

nai<T>  Receiver ( const  3ignai<T>  receive,  const  Signal<T>  chipseq,  const  long  NB,  const  long  3. con 
T  limit) 

Signal<T>  scat(3); 
stat=receive; 

//implement  receiver  for  hard  limiting 
for  (int  i=0 ;  i<s tat .lioof Columns  { )  ;  i-»-+) 

{ 

if ( { fabs (stat (0 , i ) ) )  >  limit)  stat(0,i)  =  0.0; 

} 

//  End  of  hard  limiting 

//  Implement  decision 
stat. Indmult (chipseq) ; 

stat . Reshape (NB, 3) ; 


//cout«stat; 
stat=stat .Sum( ) ; 


//Make  decisions  based  on  values  of  the  elements  of  stat 
for (i=0;i<s tat.Noof Columns ( ) ; 

{ 

if  (stat(0,i)  <=  0.0)  stat(0,i)  =  *-1.0; 
else  stat(0,i}  =  1.0; 

} 

return  stat; 
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'  Creaced  By  Abhijic  Talachi  for  ECE  449 
/  Simuiacion  of  DSS3  Communi cat ions  System 

:  f  Header  file  for  the  definition  of  the  class  Signal 
'  /  used  in  simulation  of  a  Direct  Sequence  Spread  Spectrum 
communication  system. 

//  This  header  file  defines  the  class  signal  which  is  derived  from  the 
Matrix  template  class. 

•■/Include  the  Matrix  Teplate  Definition  file 

#ifndef  SIGNAL^H 
^define  SIGNAL_H 

#inciude “Matrix. h" 

//Macro  for  uniform  random  generator 

5^define  unirand {)  (float)  (  (double)  rand () /RANB_'iAX) 

const  double  PI  =  3.1415927; 


const  int  NAMELEN  =  20;  //  Longest  possible  name 

template  <class  T>  class  Signal  ;  public  Matrix<T> 

{ 

private: 

//  A  flag  that  is  set  once  the  maximum  and  minimum  are  found 
//  This  allows  us  to  find  the  bound  only  once  if  they  are 
//  required 
short  fndbound; 

static  const  int  MU_MAX;  //  Maximum  value  of 

static  const  int  SIG_MAX;  //  Maximum  Value  of 

//  static  const  int  NAMELEN;  if  Longest  possible 

protected: 

//  The  maximum  value  of  the  signal 
T  max; 

//  The  minimum  value  of  the  signal 
T  min  ; 

//  Stores  the  name  of  the  signal 
char  tag [NAMELEN] ; 

//  Sets  the  naune  (Only  available  to  derived  classes) 
void  name  (char*  identity)  (strncpyi  tag,  identity,  NAMELEN)  ;  } 

public : 

//  Default  constructor 

SignaKconst  long  Length, const  T  mini  =  0, const  T  maxi  =  0); 

//  Converting  Matrix  to  signal 
SignaKconst  Matrix<T>  mat)  ; 

void  findbounds ( ) ;  //  Find  the  boundaries 

T  sigmax();  //  Set  the  maximum  if  necessary  and  provide  to  user 

T  sigmin();  //  Set  the  minimum  if  necessary  and  provide  to  user 

//  Returns  the  name  of  the  signal 

virtual  char*  identify ()  {return  tag;} 

//  Function  zo  change  the  amplitude  by  a  fixed  factor 
virtual  v'oid  attenuate  (const  double  factor); 

//  The  basic  types  of  Noise  models 

//  Gaussian  •  i 

void  gaussgen ( const  long  length,  const  double  mu  =  0.0,  const  double  sigma  =  1.0); 

//  Uniform 

void  unifogen ( const  long  length,  const  double  mu  =  0.0); 

i I  Swaps  first  half  and  second  half  of  signal  (similar  to  matlab) 
void  f f tshif t ( ) ; 


MU 

Sigma 

name 


50 


cacch 

i 

couc<<endi<<''Gaussian  Generator  :  Arguments  are  to  large  for  this  simulation"; 
exit iHXIT_FAILURH) ; 


double  temp=:0.0; 

double  3=sigma  /  3 . 33 ;/ /Adjusted  Sigma 

foriint  i  =  0,  n  =  12;  i  <  length;  i++,  n=i2) 

{ 

for{;  n  >  0;  n--; 

{ 

temo  +=  unirand {);/ /Call  uniform  random  generator 

} 

temp  =(3  ^  (temp  -  5)); 
data  r i ] =temD; 

} 


Generate  a  Uniformly  Distributed  signal 
template<class  T> 

void  Signal<T> : runifogen (const  long  length,  const  double  mu) 

{ 

try 

{ 

if (mu  >  MU^MAX) 

throw  "error"; 

) 

catch  (...) 

{ 

cout<<endl<< "Uniform  Generator  ;  Arguments  are  to  large  for  this  simulation"; 
exit (£XIT_FAILURE) ; 

} 

//  Generate  a  uniform  random  distribution  with  mu  as  the  mean 
for(int  1=0;  i  <  length;  i++) 

{ 

dataCi]  =  (T)(rand()  «  0.5  +  mu);//Make  it  uniform  distribution  with  meein  =  mu 

} 

} 

.'/  Swaps  the  first  half  of  the  signal  with  the  second  half 
Similar  to  the  function  in  matlab  which  does  the  sa.me. 
tempiate<class  T>void  Signal<T> : : f f tshif t  ( ) 

int  midot; 
if (coi%2) 

{ 

midpt  =  (int)  (col  /  2)  +  1;//Center  as  defined  by  the  required  function 

} 

else 

{ 

midot  =  (col  /  2); 

} 

//Divide  the  signal  into  two  parts  and  swap  each  part 
for (int  i  =  0;  i  <  midot;  i++) 

( 

temp=data (i ] ; 

daca  [  ij  =data  [i+midpt]  ; 

data  [  i’i-.midpt  1  =temp; 

\ 


template<class  T>const  int  Signal<T> :  : MU_MAX  =  10;  //  Maximum  value  of  MU 
tempiate<class  T>const  int  Signal <T> ::  SIG_MAX  =  5;  //  Maximum  Value  of  Sigma 
/const  int  Signal NAMELEN  =  20;  //  Longest  possible  name 

=  endif 
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■  Created  by  Dhiraj  Gandhi 

,  Simulation  of  DSS3  Communication  System 

■  Matrix  class 


^ifndef  MATRIX_H  //To  avoid  multiple  inclusion  of  the  same  class 
^define  MATRIX^H 

?include<iostream.h> 

? include<stdlib.h> 

»include<math - h> 

#include<assert .h> 

#include<string . h> 

#include<time . h> 

^define  EQUAL  0 
4define  N0T5QUAL  1 


,'/  Matrix  Class  Definition 
template<class  T> 
class  Matrix 
{ 

protected: 

//  Rows  -->  No  of  rows  in  the  matrix 

//  Clximns  — >  No  of  columns  in  the  matrix 

//  Data  - >  Will  point  to  starting  address  of  array 

//  Rows  and  columns  sometimes  may  greater  than  limit  of  int  hence  chosen  as  long 

long  rows , columns ; 

T*  data; 


public: 

//  Default  constructor  with  long  row  .  columns  and  value  equall  to  zero 
Matrix{const  long  row=0, const  long  col=0, const  T  value=0); 

If  Copy  constructor 
Matrix(const  Matrix<T>&  mat)  ; 

//  Constructor  taking  the  char’  as  a  arguemnet  and  converts  into  the 
//  required  type  .  ^  ^ 

//  This  constructor  can  convert  {"2  2  2,2  2  2,2  2  2  "}  into  3x3  matrix  of  2 
Matrix (const  char  ’  str) ; 

//  Destructor 
-Matrix ( ) 

{ 

delete  (]data; 

} 

//  Return  Total  No  of  Rows  of  the  matrix 
long  NoofRows ( ) 

(  return  rows;  } 

//  Returns  Total  no  of  columns  of  the  matrix 
long  NoofColumns ( ) 

(  return  columns;  } 

//  Overloading  ()  operator  for  indexing 
T  ooerator ()  (const  long  row, const  long  coDconst 
{ 

try 

if  (row>=:rows  &&  row  <0  col  >=columns  &&  coi<0) 
throw  "error"; 

} 

catch  (...) 

cout<<endl<< "  Row  no. and  Column  no  should  be  less  than  \ 
total  rows  and  columns “ <<endl ; 

return  0; 

) 

return  data [ (row’columns) +col ] ; 


52 


/  Overloading  )  operator  for  indexing  (reference) 

TSc  operator  ()  (const  long  row,  const  long  col) 

{ 

try 

{ 

if(row>=rows  ScSc  row  <0  Sci  coi  >=coluinns  &&  col<0) 
throw  ■’error*'; 

} 

catch  (...) 

C 

cout<<endl<< "  Row  no. and  Column  no  should  be  less  than  \ 
total  rows  and  columns " <<endl ; 
return  data [0] ; 


; 

return  datat (row*columns ) +coi ] ; 


//  Overloading  composite  operators 
Matrix<T>&  operacor+=  (const  Matrix<T>  rat)  ; 
Matrix<T>&  operator-=  (const  Matrix<T>  mat)  ; 


//  Overloading  composite  operators  for  value 
Matrix<T>Sc  operator +=  (const  T  value)  ; 
Matrix<T>&  operator-=  (const  T  value)  ; 
Matrix<T>&  operator '=  (const  T  value)  ; 
Matrix<T>&  operator/ =  (const  T  value)  ; 


//  Overlaoding  binary  operators 
Matrix<T>  operator- ( ) const ; 
Matrix<T>  operator+  (Matrix<T>  mat) 
Matrix<T>  operator- 
Matrix<T>  operator'* 

Matrix<T>  operator+ 

Matrix<T>  operator- 
Matrix<T>  operator* 

Matrix<T>  operator/ 


(Matrix<T>  mat) ; 

(const  Matrix<T>  mac) const, 
(const  T  value) const ; 

(const  T  value) const; 

(const  T  value) const ; 

(const  T  value) const; 


//  Overloadind  =  operator 

Matrix<T>&  opera tor =  (const  Matrix<T>  mat)  ; 
Matrix<T>Se  opera tor=  (const  T  value)  ; 


//  Overloading  comparisons  operator  for  matrix 
short  operator  ==  (const  Matrix<T>  mat) const; 
short  operator  !=  (const  Matrix<T>  mat)  const; 

//  Overloading  comparison  operators  for  value 
short  operator  ==  (const  T  value) const ; 
short  operator  !=  (const  T  value ) const ; 

//  Multiplies  individual  elements  of  the  matrix 
void  Indmult (const  Matrix<T>  mat)  ; 


//  Calculate  inverse  of  given  matrix 

Matrix<T>  Inverse{); 

//  Calculate  the  determinant  value  of  given  square  matrix 

T  Det  ( )  ; 

//  Returns  the  submatrix  part  of  gU-en  matrix  specified  by 

//  startrow,  started  and  endrow,  ended 

Matrix<T>  SubMacrix( const  long  startrow.  const  long  started,  \ 

const  long  endrow,  const  long  endcol) ; 

/  /  Reshapes  the  matrix  with  new  rows  and  columns 
void  Reshape (const  long  newrow,  const  long  newcol); 

//Returns  the  signs  of  the  each  element  of  the  Matrix 
void  Sign(); 

//  Converts  the  Matrix  into  the  row  Matrix 
void  ConvertToRowMatrix ( )  ; 
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/ /  Converts  the  Matrix  into  the  Column  Matrix 
void  ConvertTcColumnMatrixO ; 

■/  Swaps  the  Column  first  with  Column  j econd 
void  SwapColumns ( cons t  long  rirst, const  *ong  second); 

'/  Swaps  the  Row  first  with  Row  second 
void  SwapRows ( cons t  long  first, const  long  second); 

/ /  Discrete  Fourier  Transform 
Matrix<T>  DftO; 

//  Convolution  function 

Matrix<T>  Convo( const  Matrix<T>  mat); 

' !  Sum  of  the  all  elements  of  given  matrix 
j'Iatrix<T>  SumO; 

//  Finding  the  no  of  errors  between  two  given  matrices 
long  NoofErrors (const  Matrix<T>  mat) ; 

//  Polymorphic  function  for  addition 
virtual  Matrix<T>  add(const  Matrix<T>  mat)  ; 

//  Overloading  «  and  >>  operators 

friend  ostream&  operator<< (ostreamfic  os ,Matrix<T>&  cstr) ;  //  Overloading  <<  operator 
friend  istream&  operator>> (istream&  is ,2-Iatrix<T>&  cstr);  //  Overloading  »  operator 


; ; 


Default  constructor  with  long  row  ,  columns  and  value  equall  to  zero 
template<class  T>Matrix<T> :: Matrix (const  long  row, const  long  col, const  T  value) 

{ 

rows = row; 
columns=col ; 
try 
{ 

if((data=new  T [rows* columns ]) s-NULL) 
throw  "error"; 

} 

catch  (...) 

cout<<endl<< "  Memory  Allocation  Error " <<endl ; 

//  temp  is  temporary  variable  that  we  can  use  to  avoid  same  multiplications  in  the  loop 
long  temp=rows * columns ; 
for (long  i=0 ; i<temp; i++ ) 

data[i] =value; 


//  Copy  constructor 

template<class  T>Matrix<T> :: Matrix (const  Matrix<T>  &mat) 

{ 

rows  =mat . rows ; 
columns=mat . columns ; 
long  temp=rows ’columns ; 
try 
{ 

if{(data=new  T  [  temp]  )  ==NULL) 
throw  " error " ; 


catch  (...) 

{ 

cout«endl«'*  Memory  Allocation  Error"<<endl; 

} 

for (long  i=0 ; i<temp; i++) 

data  C i ] =mat . data ( i ] ; 


/  Constructor  taking  char  *  as  arguement  .  ^ 

/  This  constructor  can  convert  ("2  2  2,2  2  2,2  2  2  ")  into  3x3  matrix  of  2 
tempiate<class  T>Matrix<T> : :Matrix (const  char  *str) 

^  inc  long  MAXDIGIT=4;  //  No  of  max  digits  of  the  number 
long  length=strien (str) ; 
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long  rowcnc=0 ,  coluinncnt=0 ; 
coiuiiins=0  ; 


//  This  loop  will  determine  total  number  of  columns  and  rows 
for (long  i=b; i<length;i++) 

{ 

if(str[il==^  ') 
columncnt++; 
if  {str[il==^  ' ) 

{ 

rowcnt++; 
columncnt++; 
if {columns  ==0) { 

columns=coluinncnt ; 

} 

try 

{ 

if (columns  !=  columncnt) 
throw  *'  error  “  ; 

} 

catch ( . . . ) 

C 

cout«endl«"  No  of  columns  are  not  equal  for  every  row‘*«endl; 

} 

columncnt=0; 

} 

} 

rows=rowcnt+l ; 

data=  new  T  ( rows* columns ] ; 

long  datacnt , negative=l , convertcnt=0 ; 

T  value=0 , convert [MAXDIGIT] ; 
for(int  jpO; j<MAXDIGIT; j++) 
convert [j ] =0; 


//  This  will  take  each  individual  member  and  convert  it  into  the  required  type 
for ( i=0 , datacnt=0 ; i<length; i++ ) 

{ 

if  (strCi]==:'  '  II  str[i]  == 

{ 

value=0 ; 

for(int  m=MAXDIGIT-l ; convert [m]  ==0  &&  m  >=0;m — ); 
for ( j=0; j<=m; j++) 

{ 

value+s (convert[ j ] *pow(10, (m-j ) ) )  ; 
convert  ( j  ]  =0  ; 

} 

cout<<endl; 

da ta( datacnt 1 =negative*value ; 

datacnt++; 

convertcnt=0 ; 

negative=l; 

continue; 

if (str(i]==^-'&&{str[i-ll=='  '  ||  str[i-l]  == 

( 

negative=-l; 

continue; 

} 

convert ( conver tent ] =s  tr [ i ] - ' 0 ' ; 
convertcnt++ ; 

} 


} 

•■  /  Overloading  +=  operator  for  matrix 

cemplate<ciass  T>Matrix<T>&  Matrix<T> : : operator+=  {const  Macrix<T>  mat) 

( 

long  t emp= rows* columns ; 
for (long  i=0 ; i<temp; i++) 

data [ i ] +=mat . data [ i ] ; ; 
return  *this; 


} 
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!  Overloading  -=  operator  for  matrix 
tempiate<ciass  T>Matrix<T>&  Matrix<T>::  operator-^  (const  Matrix<T>  mat) 

long  temp=rows * columns ; 
for  (long  i=0;  i<temp;  i+-^) 

data  t i I -=mac . data [ i ] ; ; 
return  'this; 

} 

/  Overloading  +=  operator  for  value 
template<class  T>Matrix<T>&  Matrix<T> : :operator+=  (const  T  value) 

f 

long  temp^rows * columns ; 
for (long  i=0; i<temp; i++) 
data  C i ] +=value ; 
return  ^this; 


: I  Overloading  -=  operator  for  value 

template<class  T>Matrix<T>&  Matrix<T> : : operator-=  (const  T  value) 

{ 

long  temp=rows * columns ; 
for (long  i=0; i<temp; i++) 
data (i }-= value; 
return  ♦this; 


} 

//  Overloading  ♦=  operator  for  value 

template<class  T>Matrix<T>&  Matrix<T>: : operator* =  (const  T  value) 

C 

long  temp=rows*columns ; 
for (long  i=0 ; i<temp; i++) 
data(i] *=value; 
return  *this; 


) 

//  Overloading  /=  operator  for  value 

tsmplate<class  T>Matrix<T>&  Matrix<T> :: operator /=  (const  T  value) 
{ 

long  temp=rows * columns ; 
for (long  i=0; i<temp; i++) 
data[i] /=value; 
return  *this; 


} 

//  Oveloading  Unary  -  operator 

template<class  T>Matrix<T>  Matrix<T> operator- () const 

{ 

Matrix<T>  result (rows, columns) ; 
for (long  i=0; i< (rows*columns) ;i++) 
result -data [i] =-data[i] ; 
return  result; 

} 


i !  Overloading  *  operator  for  Matrix 

template<ciass  T>Matrix<T>  Matrix<T> :: operator*  (const  Matrix<T>  mat) const 
{ 

try 

( 

if(  columns  !=  mat. rows) 
throw  •*  error  “  ; 

} 

catch  (...) 

^  couc<<endl<<“  Multiplication:  Columns  of  first  Matrix  =  Rows  of  second  matrix'*«endl 
return  Matrix<T>(); 

) 

long  temp; 

Matr ix<T>  result ( rows , mac . columns ) ; 
for  (long  i=0;  i<rows;  i++) 

{ 

cemp=i ‘columns ; 
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for  (long  j=0;  j <mac. columns ;j++) 

{ 

for  (long  k=0;  ;« columns;  )c++) 

{ 

result(i.j)  ■>•=  daca( cemp-nkl  *  mat(k,j); 

} 

} 

i 

return  result; 


■  Overlading  +  operator  for  value 

cemplate<class  T>Matrix<T>  Matrix<T> :: operator+ (const  T  vaiue)consc 

Matr ix<T>  resul t ( rows , columns ) ; 
long  temp=rows*columns ; 
for (long  1=0 ; i< temp; i++) 

resul c . data [ i ] =data  l i ] +value ; 
return  result; 


•/  Overloading  -  operator  for  value 

template<class  T>Matrix<T>  Matr ix<T> operator- (const  T  value)const 

/ 

Matrix<T>  result (rows, columns)  ; 
long  temp=rows*columns ; 
for (long  i=0 ; i<temp; i++) 

resul  t .  data  [  i  ]  =sdata  [  i  ]  -  value  ; 
return  result; 

} 


/  Overloading  *  operator  for  value 
csmplate<class  T>Matrix<T>  Matrix<T> :: operator* (const  T  value) const 

Matrix<T>  result (rows, columns) ; 
long  temp=rows*columns; 
for (long  i=0; i<temp; i++) 

result . data [ i ] =data [ i ] * value ; 
return  result; 


■ /  Overloading  /  operator  for  value 

t9mplate<class  T>Matrix<T>  Matrix<T> operator/ (const  T  value) const 

Matrix<T>  result (rows, columns) ; 
long  temp=rows*columns ; 
for  (long  i=0;  i<temp;  i'»-+) 

result.data[i] =data(i] /value; 
return  result; 

} 


/  Overloading  =  operator  for  matrix 

template<class  T>Matrix<T>&  Matrix<T> : : operator=  (const  Matrix<T>  mat) 

//  If  the  size  if  the  matrix  is  not  matching  with  the  new  matrix  then 
//  delete  the  old  matrix  and  create  the  new  matrix 
if (rows  !=  mat. rows  ||  columns!=  mat . columns ) 

{ 

rows=mat .  rows  ; 
columns=mat .  columns  ; 
delete  [Idata; 
try 
{ 

if((data=  new  T [rows* columns ]) ==NULL) 
throw  “error”; 

} 

catch  (...) 
r 

w 

cout<<endl<<  ”  Memory  Allocation  error " <<endl  ; 

} 

long  temp=rows*columns 
for (long  i=0 ; i<temp; i++) 

data ( i ] =mat . data ( i  ]  ; 
return  *this; 
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throw  “error"; 


) 

catch  (...) 

{ 

cout<<endl<<"  nsolve  :  Fails  on  the  given  Matrix"; 

delete  []d; 

return  Matrix<T> ( ) ; 

} 

nrow=i-l; 

for  { j=0 ; j<rows; j++) 

( 

nrow= ( nrow+1 ) %rows ; 

result  (nrow,  i)  =*  (d+j*  (columns  +  l) +coluinns)  ; 

} 

} 

delete  []d; 
return  result; 


//  This  function  is  internally  called  by  Matrix  inverse 
//  downloaded  from  the  web 

temolate  <class  T>  long  nsolve (long  rows,T  *data) 

{ 

long  i,j,k; 
long  cols; 

T  dtemp; 

cols=rows+l; 

for  (i=0;  i<rows;  i++) 

{ 

for  { j=i ;  j<rows&&*  (data+ j  *cols+j  )  ==0 . 0 ;  j++)  ; 
if  (* (data+ j*cols+j ) ==0 . 0) 

{ 

cout«endl« “nsolve  error:  singular  matrix"; 
return  1; 

} 

if  (j!=i) 

{ 

for  (k=0;k<cols;k++) 

{ 

dtemp=* (data+i*cols+k)  ; 

*  (data+i*cols+k)  =*  (data+ j  *cols+k)  ; 

*  (data+j*cols+k)  =dteinp; 

} 

} 

for  ( j=cols-l; j>=0; j — ) 

{ 

* (data+i*cols+ j )  /=  * {data+i*cols+i) ; 

} 

for  ( j=i+l; j<rows; j++) 

{ 

for  (k=cols-l;k>=i;k--) 

* (data+ j *cols+k) -=* (data+ j *cols+i)  *  * (data+i*cols+k) 

} 

} 

for  (i=rows~2 ; i>=0 ; i — ) 

{ 

for  ( j=cols-2; j>i; j — ) 

{ 

* (data+i*cols+cols-l) -=  \ 

* (data+i*cols  + j )  *  * {data+ j *cols+cols-l)  ; 

* (data+i*cols+ j ) =0 ; 

} 

} 

return  0; 

} 


//  finding  the  value  of  determinant 
template<class  T>  T  Matrix<T>::  Det ( ) 
{ 

T  pi, p2 , p3 ,d; 
long  i,j,k; 
try 
{ 


if  (columns ! =rows ) 
throw  "error"; 


} 

catch ( . .  . ) 
{ 
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ccuc«endl<<“  Zet:  "he  Matrix  should  be  square  macrix'‘<<endl 
return  0; 


for  ( i=0  ;  i<coluinns ;  i  +  + ) 

?l=p2=l; 
p3=data[i] ; 

'<=i ; 

for  ( j=l;  j<coiuinns; 

k=  (k+1 )  %coluinns  ; 
pl*=data{ { j  ^columns) +k] ; 
p2*=data[  {columns- j)  •'columns+k]  ; 

} 

?3*=(pi-p2) r 
d-f=o3  ; 

} 

return  d; 


.  /  Returning  the  submatrix  within  the  matrix  specified  by  the  start  and  end  rows  columns 
template<class  T>  Matrix<T>  Matrix<T> :: SubMatrix {const  long  startrow,  const  long  startcol 

, const  long  endrow, const  long  endcol) 


try 

{ 

if(  endrow  <  startrow  j|  endcol  <  startcol) 
throw  " error " ; 

} 

catch  (...) 

{ 

cout<<endl« "  Submatrix:  Ending  row  and  ending  column  should  be  greater  \ 
than  stating  row  and  starting  col"«endl; 

return  Matrix<T> ( ) ; 

} 

Matrix<T>  result ( (endrow-startrow+1) , endcol -startcol+1) ; 
for ( long  i=s tar trow, p^O ; i<=endrow; i++ , p++ ) 

for (long  j^startcol , q=0 ; j<=endcol ; j++,q++) 
result {p,q) =data [ i*rows+ j ] ; 
return  result; 


//  Reshapes  the  matrix  with  new  rows  and  new  columns 

te.mpiate<class  T>void  Matrix<T> :: Reshape ( const  long  newrow,  const  long  newcol) 

try 

{ 

if ( (newrow* newcol ) !  =  ( rows  *columns ) ) 
throw  "  error  '*  ; 

} 

catch  (...) 

( 

cout«endl<<  ”  Reshape:  Total  no  of  elements  in  old  and  new  matrix  should  match ''«endl 


} 

T*  dummy  =  new  T [rows* columns }  ; 
long  temp=newrow*newcol ; 

//  This  Reshape  routine  will  store  the  elements  columnwise  as  desired  by  Matlab 
for (long  i=0, k=0;i<columns; i+-) 

fordong  j=0 ;  j<t:emp ;  j+=coiumns  ,  k++ ) 
dummy [k] =data [ i  + j ]  ; 
for (i=0, k=0;i<newrow;i++) 

fordong  j=0;  j<temp;  j+=newrow,  k++) 
data [k] =dummy [i  + j ]  ; 


rows=newrow; 
columns=newcol ; 
delete  [} dummy; 


For  each  element  of  the  matrix  Sign  gives 
1  if  the  element  is  greater  than  zero, 

0  if  it  equals  zero 
-1  if  it  is  less  t.han  zero 
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7cut<<’'  wn=*<<wn<<endl ; 

iioac  arg,c,s; 

for (long  k=0 ;  k<len;k++) 

>.’nl=k; 

for (long  i=0 ; i<len; ) 

f 

imi=i  ; 

arg=vm*kml*inil ; 
c=cos  <arg) ; 
s=sin (arg) ; 

result .data [k] +=data[i] *c; 

} 

return  result; 


•Convolving  the  given  matrix  with  provided  matrix 
cempiate<ciass  T>  Matrix<T>  Matrix<T> : :Convo (const  Matrix<T>  mat) 

long  lenl=rows*columns ; 

long  1  en2=niat.  rows  *mat.  columns,  k; 

Matrix<T>  result(l, (lenl+len2-l) ) ; 
for (long  i=0; i<lenl ; i++) 

{ 

for (long  j=0 ; j<len2 ; j++) 

{ 

k=i+j; 

r esul t . data ( k] +=data [ i ] *mat . data [ 3 ] ; 

} 

} 

return  result; 

> 


■'/  Sum  all  the  elements  of  given  matrix  Columwise 
template<class  T>  Matrix<T>  Matrix<T> : : Sum( ) 

Matrix<T>  result ( 1 , columns , 0 ) ; 
long  j  ; 

for (long  i=0 ; i<columns ; i++) 
for ( j=0; j<rows; j++) 

result.datafi] +=data[ j*columns+i] ; 
return  result; 


/  Finding  the  no  of  errors  between  two  given  matrices 
oemplate  <class  T>  long  Matrix<T> : :Noof Errors (const  Matrix<T>  mat) 


long  errors-0; 
try 

if  { (mat  .rows *mat  -  columns)  !  =  (rows* columns)  ) 
throw  -error**; 

} 

catch  (...) 

^  cout«endl«  "  Noof Errors:  Total  no  of  elements  in  old  and  new  matrix  should  match- <<end 
return  errors; 

} 

long  tempsrows * columns ; 

for (long  i=0 ; i<temp; i++) 

if{data[i3  !=  mat.data[i]) 
errors++; 
return  errors; 


••  Overloading  <<  operator 

cempiate  <ciass  T>  ostream&  operator<< (ostream&  cs , Matrix<T>&  mat) 
i=0; i<mat .rows; i++) 


for ( long 

{ 
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Created  3y  Dhiraj  Gandhi  i  Abhijii  Taiachi  for  ECE  449 
Simuiacion  of  DSSS  Coinmunicauions  Syscem 
Main  singulation  program 

= include  "bpsk.h" 

*  include  ■*  inter ferer.h" 

’include  " noise. h" 

-■  constructs  the  final  received  signal 
template<class  T> 

Signal<T>  Receiver (const  Signal<T>  receive,  const  Signal<T>  chipseq,  const  long  NB,  const  long  3,  con 
St  T  limit) 

3ignai<T>  stat(B); 

stat=receive; 

//implement  receiver  for  hard  limiting 
for (int  i=0 ; i<stat .NoofColumns ( )  ; 

{ 

if ( ( fabs (stat (0 , i ) ) )  >  limit)  stat(0,i)  =  0.0; 

} 

//  End  of  hard  limiting 

//  Implement  decision 
stat . Indmult (chips eq) ; 

stat. Reshape ( NB , B ) ; 


//cout<<stat; 
stat=stat .Sum( ) ; 


//Make  decisions  based  on  values  of  the  elements  of  stat 
f or  (i=0 ; i<s tat  .Noof Columns  ( )  ;  i-c-*-) 

C 

if  (stat(0,i)  <=  0.0)  stat(0,i)  =  -1.0; 
else  scat(0,i)  =  1.0; 

} 

return  stat; 

> 

int  main( ) 

const  int  B  =  1028;  //  Symbols  transmitted  and  received 
const  int  NC  =  4;  //  Samples  per  chip 

const  int  PG  =  16;  //  processing  gain 

int  ISR;  //Interference  to  signal  Ratio 

//  To  store  the  coded  signal 
Signal <double>  Codedsig (B*NC*PG) ; 

//  Type  signal  to  store  the  received  signal 
Signal<double>  Received (B*NC*PG) ; 

cout«"This  is  a  Demo  for  specified  SNR  and  ISR  *<<endl; 
cout«“ Errors  are  calculated  for  specified  values  of  ISR  ’•«endl; 


int  SNR; 

do 

{ 

cout<<endl<< “  Enter  the  Signal  to  noise  ratio (  0  -  100  dB)  ;  "  ; 
cin>>SNR; 

}while(SNR  <0  ||  SNR  >  100); 

do 

{ 

cout<<endi<< "  Enter  the  Interference  to  signal  ratio  (0-  100)  : 
cin>>ISR; 

}while(ISR  <0  ||  ISR  >  100); 


double  Noisefactor  =  pow( 10 , (SNR/10)); 
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///  Create  a  BPSK  signal 
3PSK<doubie>  Transsig(B); 

//  To  score  the  chip  sequence 
Signal<doubie>  chips (PG); 

//  Codedsig  is  the  signal  to  be  transmitted. 

//  It  is  created  from  Transsig  by  applying  Discrete  Spread  Spectrum  techniques 
//  ,  Chips  is  Che  chips  sequence. 

chips=:Transsig .dsss  (Codedsig, NC,  3,  PG,  1)  ; 

//  Generate  a  Continuous  Wave  Interferer 
Interf erer<double>  *Int; 

//  Generate  a  Gaussian  Noise  Sequence. 

Gaussnoise<double>  GN(B*NC*PG) ; 

//Multiply  the  Noise  by  the  Reciprocal  of  the  SNR 
GN  *=  (l/Noisefactor ) ; 

cout«’*Choose  the  type  of  interferer "<<endl ; 
cout«'*l  :  Continuous  Wave  Interferer“<<endl; 
cout«"2  :  Gaussian  Interferer " «endl ; 
cout«''3  :  Uniform  Interferer “«endl; 

short  choice  =  0; 
while (choice  >3  j|  choice  <1) 
cin»choice; 

Contwave<double>  CW(B^NC^PG,  ISR,NC*3)  ; 

Gauss inter<double>  GI{B*NC*PG,ISR,0,1); 

Unifointer<double>  UI (B*NC*PG, ISR,  0)  ; 

//  Get  the  interf ererence  and  add  the  noise 
switch  (choice) 

{ 

case  1 : 

{ 

Int=  &CW; 
break; 

} 

case  2 : 

{ 

Int  =  &GI; 
break ; 

} 

case  3 : 

{ 

Int  =  &UI; 
break; 

} 

} 

//  Add  will  polymorphically  add  the  signal  ,  noise  and  chosen  interference 
Received  = (Int ->add( Codedsig) ). add (GN) ; 

//  Final  will  store  the  an  estimation  of  the  original  signal  at  the  receiver 
Signal<double>  final (B); 

//  Make  decisions  based  on  the  signal  and  thus  attain  the  final  estimated  signal 
final=Receiver (Received, chips , 64, B, 1000.0)  ; 

//  Calculate  the  difference  between  the  transmitted  and  received  signal 
long  cnt=Transsig. Noof Errors ( final)  ; 

cout«"The  result  for  SNR  ”«SNR«"db(s)  is  “«cnc«'  Errors “«endl;/ /Result  of  the  simulation 
return  0; 

} 
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